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Abstract 

We consider unitary random matrix ensembles Z~\ t e~ ntT v "^ M ^dM on the space of 
Hermitian n x n matrices M, where the confining potential V a> t is such that the limiting 
mean density of eigenvalues (as n — > oo and s,t — > 0) vanishes like a power 5/2 at a 
(singular) endpoint of its support. The main purpose of this paper is to prove universality 
CN| ' of the eigenvalue correlation kernel in a double scaling limit. The limiting kernel is built 

out of functions associated with a special solution of the Pf equation, which is a fourth 
order analogue of the Painleve I equation. In order to prove our result, we use the well- 
known connection between the eigenvalue correlation kernel and the Riemann-Hilbert (RH) 
problem for orthogonal polynomials, together with the Deift/Zhou steepest descent method 
to analyze the RH problem asymptotically. The key step in the asymptotic analysis will be 
the construction of a parametrix near the singular endpoint, for which we use the model RH 
problem for the special solution of the Pf equation. 

In addition, the RH method allows us to determine the asymptotics (in a double scaling 
. limit) of the recurrence coefficients of the orthogonal polynomials with respect to the varying 

weights e~ nVa '* on R. The special solution of the Pf equation pops up in the n~ 2 / 7 -term of 
I . the asymptotics. 

1 Introduction and statement of results 
1.1 Unitary random matrix ensembles 

On the space TC n of Hermitian nxn matrices M, we consider for n £ N and s,t£l the unitary 
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random matrix ensemble, 

J_ e -ntrV s , t (M) dM _ ( L -q 
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n,s,t 



Here, Z ns % is a normalization constant and the confining potential V S} t is a real analytic function, 
depending on two parameters s,t E R, satisfying the asymptotic condition, 

lim - — r = +00, uniformly for s,t € [— do, do] for some So > 0. (1-2) 

x'^ioo \og(x z + 1) 

Then, 

Z nst = I e~ ntTV ^ M UM 
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is convergent as n — > 00 so that the random matrix model is well-defined. 

It is well-known, see e.g. 25 , that an important role in the study of the unitary random 
matrix ensemble is played by the following scalar 2-point (correlation) kernel, 
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constructed out of the orthonormal polynomials 

with respect to the varying weights e~ nVa - t on R. Indeed, the correlations between the eigenvalues 
of M can be written in terms of the correlation kernel. More precisely, the m-point correlation 
function TZn'm satisfies (23] , 

K<g${x u ...,x m ) = det (K^ix^xj)) . (1.4) 

Further, the limiting mean eigenvalue distribution fj, Sj t has a density p s ,t which can be retrieved 
from the correlation kernel as follows, 

p s>t (x) = lim -K^\x,x). (1.5) 

n— >oo n 

The limiting mean eigenvalue distribution p, St t equals ^U] the equilibrium measure in external 
field V S f This is the unique measure minimizing the logarithmic energy 



WsM = ff log j^-^dv(x)dv(y) + J V s ,t(y)dfi(y), (1.6) 

among all probability measures fx on R. Furthermore, there exists a real analytic function q s t, 
such that 0, 

7T V ' 

where g^t denotes the negative part of i.e. = gj~ t — q~ t , with > and q^ t q^ t = 0. 
Due to condition (|1.2j) we have that q s ,t(x) — > +oo as i -> ±oo, so that fi s ^ is supported on a 
finite union of intervals, which we denote by S s f It is known |2H] that the equilibrium measure 
fj, s t satisfies the following Euler-Lagrange variational conditions: there exists a constant k s j € R 
such that 

2 J log |x - u\dns,t{u) - V S)t (x) = K S j, for x £ S s>t , (1.8) 

2 y log |x - u\dfx Stt (u) - V Stt (x) < K Sjt , for x G R \ S S)t . (1.9) 

The external field V^^ is called regular if strict inequality in (|1.9|1 holds, if the density p s> t 
does not vanish in the interior of the support § S) t, and if q Sl t has a simple zero at each of the 
endpoints of the support S St t- If one of these conditions is not valid, V S; t is called singular. The 
singular points x* are classified as follows, see |lf)[ 12 lj : 



(i) x* G R \ S s> t is a type I singular point if equality in ()1.9|) holds. Then, x* is a zero of qf t 
of multiplicity 4m with m G N. 

(ii) x* G S S) t is a type II singular point if it is an interior point of S S) t where the equilibrium 
density p S) t vanishes. Then, x* is a zero of q~ t of multiplicity 4m. 

(iii) x* is a type III singular point if it is an endpoint of the support S S) t and a zero of q Sj t of 
multiplicity larger than one. Then, x* is a zero of q s t of multiplicity 4m + 1, which means 
that p Stt (x) ~ c|x - x *|(4m+l)/2_ 



2 



In this paper, we consider external fields V s j which are such that in the critical case s = t = 0, 
= ^0,0 has a type III singular (edge) point x* with m = 1, i.e. 

Po,o(x) ~ c\x - x*| 5/2 , as x -> x*. (1-10) 

Further, we take V s> t of the special form, 

V s , t = V + sV 1 + tV 2 , (1.11) 

where Vi is an arbitrary real analytic function, while V2 is real analytic and in addition satisfies 
some critical condition which we will specify in Section fl .41 below. 

1.2 Universality in random matrix theory 

Consider for now unitary random matrix ensembles Z~ 1 e~ ntrV ^ M ^dM on the space of Hermitian 
n x n matrices M. Scaling limits of the associated correlation kernel K n show universal behavior. 

Near regular points, universality results have been established in P3 02 ll()| lll| I27j . For 
example, if x* lies in the bulk of the spectrum (i.e. x* is such that it lies in the interior of the 
support § of the equilibrium measure in external field V, and such that the equilibrium density 
p does not vanish at x*) there is a constant c such that 

1 v / * u * v \ sin7r(n — v) „„. 
lim — K n [x* + — ,x* + — )= p-i r-i. (1.12) 

n—>oo cn V cn cn/ 7r(u — t>) 

On the other hand, if x* is a regular edge point of the spectrum (i.e. x* is an endpoint of § and 
p vanishes like a square root at x*), there is a constant c such that 

lim 1 K „ ( x . + AiMAi'M-AiWAi'M 13) 

n^oo cn z / J V cn ' J cnrl a J u — v 

where Ai is the Airy function. 

Near singular points, similar results hold. In those singular cases it is interesting to consider 
double scaling limits where the external field V depends on additional parameters. In [21EJIBJI2H]) 
an external field V was considered such that there is a type II singular (interior) point x* with 
m = 1, i.e. 

p(x) ~ c(x — x*) 2 , as x — > x*. 

If an additional parameter is included in the external field, Vt = V/t, one observes for t close to 
1 the transition where two intervals in the support of the limiting mean density of eigenvalues 
merge to one interval through the critical case of a type II singular point. In the double scaling 
limit where n — > 00 and t — > 1 in such a way that c n 2 / 3 (i - 1) s G M for some appropriately 
chosen constant Co, there exists a constant c such that (for the associated correlation kernel 

K n ,t), 

Here, if crit ' n (u , u; s) is built out of functions associated with the Hastings-McLeod solution ^B] 
of the second Painleve equation. 

The main purpose of this paper is to obtain, for the random matrix models in Section 11.11 
above, a similar result near the type III singular (edge) point of Vq with m = 1. We take a 
double scaling limit (n — * 00 and s,t — * 0), and the limiting kernel if cnt > ni will be built out 
of functions which are associated with a special solution of the fourth order analogue of the 
Painleve I equation. The case of a type III singular (edge) point was also studied in the Physics 
literature [3 15 • 

In addition, the techniques that we use to prove this allow us to determine the asymptotics 
(in a double scaling limit) of the recurrence coefficients in the three-term recurrence relation 
satisfied by the orthogonal polynomials p^ ,s ^ with respect to the varying weights e -nl/s '* on M. 
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1.3 ^-functions associated with a special solution of the Pf equation 

We consider the following differential equation for y = y(s,t), which we denote as the Pf 
equation, 

s = ty- Qy 3 + ^{y 2 s + 2yy ss ) + ^Vssss^j • (1.14) 

For t = 0, this equation is the second member in the Painleve I hierarchy |20| I23j . The Pf 
equation has been studied for example in (for t = 0) and [TIE] (for general t). The 

Lax pair for the Pf equation is the linear system of differential equations 

- = U^ - = WV, (1.15) 
-4y s ( - (12yy s + y sss ) 8( 2 + 8y( + {12y 2 + 2y ss - 120i)\ 




U 2 i % S C + (12yy s + y s 



;i.l6) 



Usx = 8C 3 - 8yC 2 - (% 2 + 2y ss + 120t)C + (16?/ 3 - 2y 2 s + Ayy ss + 240s), (1.17) 

and 

W=(, V IV ( L18 ) 



C-2y 

The system of differential equations Hl.15JI - Ql.18jl can only be solvable if y = y(s, t) is a solution 
to the Pf equation Q1.14J) . For different solutions y, we have different Lax pairs. 

We are interested in the special solution y which was studied in [7J Ej. This solution 
y = y(s,t) is characterized by the vanishing of its Stokes multipliers si, S2, S5, and sq, see ^1 
for details. It was shown in [7, that y has no poles for real s and t, and that it has, for fixed 
t£l, the following asymptotic behavior, 

y (s,t) = ^(6|s|) 1/3 ^-6 2 / 3 t|sr 1 / 3 + 0(|sr 1 ), ass -► ±00. (1.19) 
3 

It has been shown in |26t Appendix A] that for t = 0, y is uniquely determined by realness and 
asymptotic condition Q1.19J) . For general t we are not aware of a similar result although it is 
supported by a conjecture of Dubrovin |15j that this should hold for general t. For s,t £M, the 
Lax pair Q1.15|) - Q1.18|) associated with this special choice of y has a unique solution ( 1 ) for 
which the following limit holds, see [71 119j. 

i ff 3 f$i(C;MA e e(c, s ,t) _^ J_ / 1 \ 
c ^ 2 (C;MV v^V-V ' 

as C -» 00 with < Arg C < 6vr/7, (1.20) 
where C3 = ( -1 ) denotes the third Pauli-matrix, and where 9 is given by 

0( C ; S ,t) = -i-C 7/2 -k 3/2 + < 1/2 . (1.21) 
105 S 

The functions <3?i and $2 w ih appear below in the universal limiting correlation kernel near type 
III singular (edge) points of Vq with m = 1. 
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1.4 Statement of results 

We work under the following assumptions. 
Assumptions 1.1 

(i) We consider external fields V S} t of the form 

V s , t = V + sV l + tV 2 , (1.22) 

where Vb, V±, and V 2 are real analytic and are such that there exists a 5o > such that 
the following holds 

lim - — -^-r = +00, uniformly for s, t G [— 5q, do]. (1-23) 

|a;|->oc log(ar + 1) 

(ii) Vo is such that the equilibrium measure v§ in external field Vo is supported on one single 
interval [a, b] C M, and b is a type III singular (edge) point of Vo with m = 1. Then, vo is 
of the form [2], 

cMz) = —h (x)^/ (b -x)(x- a) X[ a ,b](%)dx, (1.24) 

with X[a,b] the indicator function of the set [a, b], and with ho real analytic and satisfying, 

h (b) = h' (b) = 0, and h' \b) > 0. (1.25) 

Furthermore, we assume that Vo has no other singular points besides b. In particular, a is 
a regular (edge) point and we then have that 

h (a) > 0. (1.26) 

(iii) V2 is such that it satisfies the critical condition 



^-—Vj(u)du = 0. (1.27) 
b — u 

Throughout the rest of this paper we let V be the neighborhood of the real line where 
Vb, Vi, V2, and ho are analytic. 

Example 1.2 The assumptions above are valid for the particular example where Vb, Vi, and V2 
are given by, 

Vo(x) = rU 4 - ^-x 3 + \x 2 + fx, Vi(x)=x, V 2 (x) = x 3 -6x. (1.28) 
20 15 5 5 

Then, the equilibrium measure z^o is supported on the interval [—2, 2] and given by 

dU ° {x) = ^( x + 2 ) 1/2 ( x - 2 ) 5/2 X[-2,2](x)dx. (1.29) 

It should be noted that a type III singular (edge) point cannot occur when Vb is a polynomial 
of degree lower than 4. 
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Example 1.3 In the continuum limit of the Toda lattice ^2]; an external field of the form 

V h:t2 (x) = (l + ti)(V (x) + t 2 x) 

has been studied. This deformation of Vo can be written in the form (jl.22|) (so that it is 
included in the class of external fields studied in this paper). Indeed, if we let V\{x) = x and 
V 2 {x) = Vq(x) + cx, with c some constant chosen such that the critical condition (|1.27|) holds, 
then 

V tl ,t 2 =Vo + 8V 1 +tV 2 , 

with s = t 2 + t±t2 — cti and t = t\. 

Remark 1.4 In Section [21 we will show that assumption (iii) is equivalent to the vanishing of 
the equilibrium density du ^^ at the right endpoint b, where v 2 is the unique measure which 
minimizes Iy 2 {y), see ()1.6|) . among all signed measures u, supported on [a,b] and having zero 
mass, v([a, b]) = 0. 

Remark 1.5 The case where the left (instead of the right) endpoint of the support is singular 
can be transformed to our case by considering the external field V s ,t(— x). 

Remark 1.6 Without giving any mathematical details, we now describe the transitions that 
can occur for s and t near 0. First, if we let t = and s vary around 0, one typically observes 
the transition from the regular one-interval case to the singular case and back to the regular 
one-interval case. Next, for s = and t around 0, we can observe the transition from the regular 
one-interval case to the regular two-interval case. Finally, letting both s and t vary around 0, 
we can observe one of the above described transitions, or the critical transition where a type II 
singular point moves to the endpoint b, where it becomes a type III singular point before moving 
on as a type I singular point. 

Further, to describe our results, we have to introduce constants c, c\, and c 2 , 

c^W^) 2/7 >0, c, =- t75 AL^L_ > (1.30) 

where ho is the real analytic function appearing in (|1.24[) . and where the functions hi and hi 
are defined as, 



hj(x) = f s/(b-u)(u-a)V-(u) , for x G [a, b] and j = 1, 2. (1.31) 

7T J a 11 X 

1.4.1 Universality of the double scaling limit 

Our main result is the following. 

Theorem 1.7 Let V s> t = Vq + sVj. + tV-2 be such that Assumvtions II. ll above are satisfied. We 
take a double scaling limit where we let n — > oo and at the same time s,t — > 0, in such a way 
that lim n 6 ^ 7 s and lim n 4//7 i exists, and put 

s = ci • limn 6 / 7 s G M, t = c 2 • limn 4/7 t G M, (1.32) 

where the constants c\ and c 2 are defined by (|1.3Uj) . Then, the 2-point kernel Kn satisfies the 
following universality result, 

lim _L ( b+ ^^ + v-s =K c*t,m^ v . SQjto)) (133) 
cn z i ' \ cn z ' ' cn z ' ' / 
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uniformly for u, v in compact subsets of R. Here, K mt,lu is built out of the functions _>x an d 
$2 defined in Section lOl 

K [u,v,s,t) = — — . (1-34) 

2m[u — v) 

Remark 1.8 Since y(s,t) has no poles [7j for s,t £ 1, the kernel K cnt,lll (u, v; s, t) exists for 
all real u,v,s, and t. Furthermore, using a similar argument as in Lemma 2.3 (ii)], one can 
show that e 7 ™/ 4 ^! and e 7 ™/ 4 -^ are real. It then follows that K CTlt ' lll (u, v; s, t) is real for real u, 
v, s, and t. 

Remark 1.9 It is possible to give an integral formula for J£ cnt,m . Using the fact that 
satisfies the second differential equation of the Lax pair (|1,15|) , we have that 

-^(C;M) = *2(C;M), and -^(C;*,*) = (C- M*, t))*i(C;M)- 

Using ()1.34j) this yields, 

g^crit.III ^ 

— — (u,v;s,t) = -— $i(u;s,i)*i(v;s,i). 

Now, since Iiin._H._oo J fir orit ' in (i*, «; a, t) = 0, which can be shown using a Deift/Zhou steepest 
descent method argument ^2!) it then follows that K CTlt ' has the following integral formula, 



K crit ' m (u, v; s, t) = — — [ $i(_; a, i) a, t)da. (1.35) 

27TZ ,/__ 

Remark 1.10 Theorem 11.71 can be generalized to the case where the support of v$ (the equili- 
brium measure in external field Vq) consists of more than one interval. Then, the proof becomes 
much more technical, although the main ideas remain the same. We comment in Remark 13. 
on the modifications that have to be made in the multi-interval case. 



1.4.2 Recurrence coefficients for orthogonal polynomials 

It is well-known |3U| that the orthonormal polynomials pk = p£"' s,t ' satisfy a three-term recur- 
rence relation of the form, 

xp k (x) = a k+1 pk+i(x) + b k p k (x) + a fc p fc _i(x), (1.36) 

where a& = a^' s '^ > and bu = &^ n,s '^ G R (we suppress the s and t dependence for brevity). 
In the generic case where Vq has no singular points, the recurrence coefficients for s = t = 
have the following asymptotics, see e.g. j__j, 

a M,o) = _Z_ + ( -( n -l) j 6 (»,o I o)__+_ + -( n -i) j asn-^oo. (1.37) 

For singular potentials Vq, the constant terms in the expansions Q1.37|) remain the same, but the 
error terms behave differently [21 _j- In our case of interest, where we have a type III singular 
(edge) point of Vq with m = 1, the error term is of order 0(n -2 / 7 ), and the coefficient of the 
n~ 2 / 7 term is expressed in terms of the special solution y of the Pf equation discussed in Section 
Ol 
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Theorem 1.11 Let V 8 % be such that Assumptions ^^ above are satisfied 
term recurrence relation ()1.36j) satisfied by the orthonormal polynomials p^ 
to the weight function e _nt/s ' ( . Then, in the double scaling limit where n - 
such a way that lim n 6 / 7 s and lim n 4 / 7 i exists, and put 

s = ci • limn 6/7 s G R, t = c 2 • limn 4/7 i G M, (1.38) 

with c\ and c 2 given by (|1.30() . we /iave 

^ = ^ + £ 2/(cm 6 / 7 s, c 2 n 4 / 7 t)n- 2 / 7 + 0(n-V\ (1.39) 
= ^ + ^»(-o, to)n- 2 / 7 (l + o(l)), 

and 

b (n,s,t) = b +^L + l y( Cl n 6 / 7 s, c 2 n i ' 7 t)n~ 2 / 1 + 0(n- 3 / 7 ) (1.40) 

= ^ + 1 ^o,to)n- 2 / 7 (l + o(l)), 
2 c 

where the constant c is given by 1)1.30)) . and where y is the special solution of the Pj equation 
discussed in Section ll~Bl 

Remark 1.12 Note that the expansions of the recurrence coefficients are of the same form 
as the conjectured (by Dubrovin ^51 Main Conjecture, Part 3], see also ^1]) expansions for 
solutions of perturbed hyperbolic equations. Here, the perturbation parameter e plays the role 
of 1/n in our context. 

Remark 1.13 For polynomials which are orthogonal on certain complex contours, it can occur 
that the equilibrium density vanishes like a power 3/2. Asymptotics of the recurrence coefficients 
in this case were obtained in ^Hj- Here, a special solution of the Painleve I equation occurs 
instead of a solution of the Pf equation and the asymptotics are in powers of n -1 / 5 . 

Observe further that in there is no term of order n -1 / 5 in the asymptotics. In 1)1.39)1 
and l)1.4U)) we see that there is no term of order n -1 / 7 . In the proof of Theorem 11.111 this term 
will drop out in a similar way as the n _1//5 -term in |16j . 



. Consider the three- 
-> oo and s, t — * 0, in 



1.5 Outline of the rest of the paper 

We prove our results by characterizing the orthogonal polynomials via the well-known 2x2 
matrix valued Fokas-Its-Kitaev Riemann-Hilbert (RH) problem |17j and applying the Deift/Zhou 
steepest descent method ^SJ to analyze this RH problem asymptotically. This approach has 
been used many times before, see e.g. pi5)l8llTU ) lrT ) IT7 ^ l ^ l3T ) l32]. 

An important step in the Deift/Zhou steepest descent method is the construction of so- 
called g- functions associated with equilibrium measures. Those equilibrium measures will be 
constructed in Section [21 In order to deal with the deformations V s> t of Vo, we use modified 
equilibrium problems where we allow the measures to be negative, which was also done in 
[3 El DEI • Another modification of the equilibrium problem is that we choose the support of the 
equilibrium measure fixed, instead of allowing it to choose its own support. 

In Section |21 we perform the Deift/Zhou steepest descent analysis to the RH problem Y for 
orthogonal polynomials. Via a series of transformations y^T^S'^i^we want to arrive 
at a RH problem for R which is normalized at infinity (i.e. R(z) — > I as z — > oo) and with 
jumps uniformly close to the identity matrix. Then, R itself is close to the identity matrix. 
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By unfolding the series of transformations we then get the asymptotics of Y. The key step in 
this method will be the local analysis near the endpoints a and b. Near the regular endpoint 
a, we construct (in Section l3.5j) a parametrix built out of Airy functions. Due to the modified 
equilibrium measures, which have a fixed support, we also need to make a technical modification 
in the construction of the Airy parametrix, compared with the parametrix as used e.g. in [S]. To 
construct the local parametrix near the singular endpoint b (in Section f3.6j) we use a model RH 
problem associated with the special solution y of the Pf equation as discussed in Section 11.31 

The results of Section |3] will be used in Section 0] to prove the universality result for the 
correlation kernel (see Theorem ll.7|) and in Section |S] to determine the asymptotics of the 
recurrence coefficients (see Theorem II .11(1 . 



2 Equilibrium measures 

We consider external fields V S) t = Vq + sV\ +tV2 which satisfy Assumptions II . ll in the beginning 
of Section fl. 41 In order to perform the Deift/Zhou steepest descent analysis to the RH problem 
for orthogonal polynomials one would expect to use the equilibrium measure // s t in external 
field Vs t minimizing iy st (/^), see p.6|) . among all probability measures [i on R. However, as in 
[HI El QH] it will be more convenient to use modified equilibrium measures v s> t which we allow 
to be negative. Furthermore, unlike in 0IB1E]) we take the support of the measures v s ^ to be 
fixed instead of letting it depend on s and t. 

The aim of this section is to find measures i/ s $ (depending on the parameters s,t G R) 
supported on the interval [a, b] C R (where [a, b] is the support of the equilibrium measure Vq in 
external field Vo), such that v Si t([a,b\) = 1, and such that they satisfy the following condition: 
there exist £ Si t G R such that for every 5 > there are e, k > sufficiently small such that for 
s,t G f-e,e|, 



(2.1) 



2 J log \x - u\dv S) t{u) - V Sjt (x) = £ Sj t, for x G [a, b]. 

2 J log \x - u\du S)t (u) - V s>t (x) < i Sjt - k, for x G R \ [a - S, b + 6]. (2.2) 

We seek v s % in the following form, 

v s ,t = + sv\ + tu 2 , (2.3) 

where vq is the equilibrium measure in external field V$ minimizing Iy (u), see ()1.6|) . among 
all probability measures v on R. From Assumption 11.11 (ii) we know that vq can be written as 
follows 

dv Q {x) = i>Q t +(x)x\a,b\{x)dx, (2.4) 

where X\a,b] i s the indicator function of the set [a, b], and where ^0,+ is the +boundary value of 
the function 

^o(z) = ^-R(z)h (z), for z G V \ [a, b], (2.5) 
with ho analytic in the neighborhood V of the real line, and with 

R(z) = ((z-a)(z-b)) 1/2 , for z G C\ [a,b]. (2.6) 

Here, we take the principal branch of the square root so that R is analytic in C \ [a, b\. Further, 
since a is a regular (edge) point and since b is a type III singular (edge) point with m = 1, we 
have, cf. (|T23|) and <Pfll . 

h (a) > 0, ho(b) = h' (b) = 0, and h'^b) > 0. (2.7) 
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Since Vq is assumed to have no other singular points besides b, we know (cf. (|1.8|) and (|1.9|> ) 
that vq satisfies the following condition: there exists ^£8 such that 

2 J log \x - u\dv a (u) - V (x) = £ , for x £ [a, b], (2.8) 

2 / log \x - u\dv (u) - Vq{x) < £ , for x £ R \ [a, b}. (2.9) 



We will now construct the two measures v\ and v%. In order to do this we introduce the 
following auxiliary (analytic) functions, 

kj( - z) = ^Tif 1 for z G V and j = X ' 2 ' (2 ' 10) 

where 7 is a positively oriented contour in V with [a, b] and z in its interior, and where R is 
given by ()2.6|) . Observe that, using the fractional residue theorem, one has, 

hJx) =-—-[ R+(u)VUu)-^-, for x e [a,b], (2.11) 

TTl J a J U-X 

where the integral is a Cauchy principal value integral. So, hj is real on [a,b]. Observe that by 
Assumption ll.il (iii) and (|2.11[) . 

h 2 {b)=0. (2.12) 

Lemma 2.1 Define two signed measures v\ and vi supported on [a,b] as 

dvjix) = i/) j>+ (x)x[ a ,b]dx, 3 = h 2 , ( 2 -!3) 

where X[a,b] * s ^ e indicator function of the set [a,b], and where is the ^boundary value of 
the function 



1 hj(z) 
2ri R(z) ' 



^( z ) = — -^7IV> M^e V\[a,6]. (2.14) 



Here, hj is given by (|2.10j) . see a/so (j2.11l) /or ite expression on [a,b], and R is given by (|2.6|) . 
Then, Vj has zero mass, i.e. 

Uj ([a,b})= j ip jt+ (u)du = 0, (2.15) 
and i/iere exist constants £j £ M. such that 

2 J log |x - n|d^(u) - Vj(x) = £j, for x G [a, 6]. (2.16) 

Proof. Define, for j = 1,2, the auxiliary functions 

1 f b du 
F ' W = 5SSW 7. ^W^'C'' — ' fa = eC \M. (2.17) 

which, by standard techniques and by (|2.1U|) and (|2.14j) . are equal to 

d£ 



\v;{z)-Tti^{z), for z£V\[a,b], 
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where 7 is a positively oriented contour in V with [a, b] and z in its interior. This, together with 
the fact that + = —^j- on (a, 6), yields 

F jj+ (x) - Fj-(x) = -2mi/)j j+ (x), for x G [a, b], (2.18) 

+ Fj-(x) = V-(x), for x G [a, b]. (2.19) 

Since Fj is analytic in C \ [a, 6] and since, by (|2.17|) . i^-(^) = 0(z -2 ) as z — ► 00, a standard 
complex analysis argument, shows that 



— [ Fj ' +{U) Fj ^ u) ds = Fj(z), for z£C \ [a, 6]. 



2m 

By (|2~TH|) . this yields, 



Fj(z) = — — du = — z / ipj )+ (u)du + 0(z ), as z — > 00. 



(i 



Comparing this with the fact that Fj(z) = C(z -2 ) as z — » 00, we obtain f ij)j t +(u)du = 0, so 
that (|2.15|) is proven. 

It remains to prove (|2.16|) . It is straightforward to check that, 

FAz) = - [ b ^' +( ^ dn = -ttMz) + -<£ t^-dt, for z G V \ [a, b], 
J,, " "- 2 / 7 £ - z 

so that, using the fractional residue theorem, 

/•& ^! ■ _|_(li) 

^',±0*0 = —Kiipj^ix) — + — du, for x G [a, b]. 

J a U — X 

From (|2.19|) and the fact that + + ipj- = on [a, b] this yields, 

^- ^2 j log \x - u\dvj(u) + Vj(x)J = 2 £ ~~ du + 0*4,+ (&) + F i -0*0) = °- ( 2 - 20 ) 

This proves (|2TTS|) . □ 
Corollary 2.2 Lei za^ = vq + s^i + if*. Then, du S)t (x) = if) s j,+(x)x\ a ,b]dx, where 

il>s,t = ik + 8ifo.+tifo, onV\[a,b], (2.21) 

with ipQ given by (|2.5jl and and ip2 given by (|2.14|) . So, za^j is supported on [a, b] and has mass 
one, i.e. v s t([a, 6]) = 1. Further, there exist constants £ s t G M such that for any 5 > i/iere are 
e, k > sufficiently small such that for s,t G [— e, e] i/te conditions (|2.1|) and (|2.2[1 are satisfied. 

Proof. Since i/ s ^ = z^o + s^i + if*, from (|2.15|) . and from the fact that vo{[a, b\) = 1 it is clear 
that v s j([a, b]) = 1. Next, with £ Sjt = £0 + s^i + U2, we have 

2 J log |x - u|di/ a , t («) - V 8 ,t(x) - e s<t = I (x) + sh(x) + t/ 2 (a?) (2.22) 

where 

Ij(x) = 2 / log \x — u\dvj(u) — Vj(x) — £j, j = 1, 2, 3. 
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Then, condition (f!TT|) follows from and Q2.lt) |) . Now, by using fiHfy and the fact that 

Io(x) — > — oo as |x| — > oo, there exists k > such that 

I < -gK, onM\[a-<y,6 + «J]. (2.23) 

Further, one can check that I\ and I2 are bounded on R \ [a — S, b + 6], and thus there exists 
e > such that for s, t G [— e, e], 

a/i + tl 2 < -k, on R \ [a — 5,b + 5]. (2.24) 

Inserting (f2~2"3)) and (|2~2^|) into (|!T22"|) we obtain condition (|^2l . □ 

Remark 2.3 The measure 1/1 (1/2) is the equilibrium measure that minimizes Iviiy) (-fyfeO")) 
among all signed measures u, supported on [a, b] with u([a, b]) = 0. The measures u s ^ on the 
other hand minimize Iv at { u ) among all signed measures supported on [a, b] with z^([a, b\) = 1. 

Observe that since has a strictly positive density on (a, b) (since uq has no type II singular 
points) we have for any 5 > that v s j is positive on (a + 5, b — 5) for s, t sufficiently small. 



3 Riemann-Hilbert analysis 

3.1 RH problem for orthogonal polynomials 

For each fixed n, s, and t, we consider the Fokas-Its-Kitaev Riemann-Hilbert problem ^7] cha- 
racterizing the orthogonal polynomials p k n ' s '^ with respect to the weight functions e~ nVs ' t . We 
seek a 2 x 2 matrix-valued function Y (z) = Y(z; n, s, t) (we suppress the n, s, and t dependence 
for brevity) that satisfies the following conditions. 



RH problem for Y: 

(a) Y : C \ R -> C 2x2 is analytic. 

(b) Y possesses continuous boundary values for x G R denoted by Y + {x) and Y_(a;), where 
l + (x) and Y"_(x) denote the limiting values of Y(z') as 2/ approaches x from above and 
below, respectively, and 

y+(x) = y_(x) ( * l , fonei. (3.1) 

(c) Y has the following asymptotic behavior at infinity 

Y{z) = {l + 0{z~ 1 )) M, asz-oo. (3.2) 



The unique solution of the RH problem is given by 

Y(z) = 

^-27riK n _ip n _i(z) -Kn-i 



- du 



u — z 
p n -i{u)e- nV ^ 



du 



for z G C \ R, (3.3) 



u — z 



where p& = p) ' is the k-th degree orthonormal polynomial with respect to the varying weight 

e" nl/s ' ! , and where Kk = k^' 5 '*'* > is the leading coefficient of pk- The solution 1)3. 3|) is due to 
Fokas, Its, and Kitaev |T7], see also |51 ITUl ITT] . 
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It is now possible to write the 2-point kernel Kn' l \ see (|1.3|) . in terms of Y. Indeed using 
the Christoffel-Darboux formula for orthogonal polynomials and the fact that det Y = 1 (which 
follows easily from 1(3.1(1 . (|3.2|) . and Liouville's theorem), we get 

K^(x,y) = e -f Vi.tWe-f %.*(«) _J_ ( i) Y±\y)Y ± (x) Q . (3.4) 

So, in order to prove Theorem 11.71 we need to analyze the RH problem for Y asymptotically. 
We do this by applying the Deift/Zhou steepest descent method [I3J to this RH problem. 

3.2 Normalization of the RH problem at infinity: Y \— > T 

In order to normalize the RH problem for Y at infinity, the equilibrium measures v s t, introduced 
in Section |21 play a key role. Consider the log-transform g s j of i/ s t, 

g s ,t(z) = / log(z - u)du s> t(u), for z G C\ (-00,6]. (3.5) 

J a 

Here, we take the principal branch of the logarithm so that g Si t is analytic in C \ (—00, b]. We 
now give properties of g s j which are crucial in the following. From H3.5j) and condition (j2.1|) it 
follows that 

g s ,t,+(x) + g s ,t ,-{x) - V s>t (x) - £ s ,t = 0, for a; € [a, 6]. (3.6) 
Another crucial property is that 

g s ,t,+(x) - g s ,t-(x) = 2-iri / du s>t (u), for (3.7) 

J X 

so that since u S) t is supported on [a,b] and has mass one (see Corollary 12.2(1 . 

I 2ttz for x <C a 

9s,t,+( x ) -9s,t,-(x) = S n ' ' (3.8) 

I 0, for x > b. 

Now, we are ready to perform the first transformation Y 1— ► T. Define the matrix valued 
function T as 

T(z) = e- 1 2 ni ^Y(z)e- n9sAz) ' 73 e \ nt ^ , for * G C \ R, (3.9) 

where £ Sj t is the constant that appears in the variational conditions ((2.1() and ((2.2(1 . and where 
03 = ( -1 ) denotes the third Pauli- matrix. Using ()3,fiJ) . (|3.S(1 . the RH conditions for Y, and the 
fact that g s ,t{z) = log 2 + 0(1/ z) as 2 — ► 00, it is straightforward to check that T is a solution 
to the following RH problem. 

RH problem for T: 

(a) T : C \ R -► C 2x2 is analytic. 

(b) T + (x) = T_(x)vt(^c) for x G R, with 

^ e -n{g a ,t,+-gs,t-) 1 \ 

, , , on (a, 6), 

r, ^ (3.10) 

'j g"(Ss,t,++9s,t,--l / s ,t-£ Sj t)\ 

Q 1 I, onR\(a,6). 
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(c) T(z) = I + 0(l/z), asz^oo. 

Remark 3.1 From l|3.7j) we see that the diagonal entries of vt on (a, b) are rapidly oscillating 
for large n. Further, using condition (\2.2\i and ()3.5|) . we see that vt — I decays exponentially on 
R\ [a-6,b + 5\. 

3.3 Opening of the lens: T y-> S 

Here, we will transform the oscillatory diagonal entries of the jump matrix vt on (a, b) into 
exponentially decaying off-diagonal entries. This step is referred to as the opening of the lens. 

Introduce a scalar function cf) Sj t as, 

</> 8tt {z) = -7T» f ip s ,t(Od£, for z G V \ (-00, 6], (3.11) 



where the path of integration does not cross the real line, and where ip s ,t is defined by (|2.21j) . 

The important feature of the function 4>s,t is that by (|3.7|) . s ,t,+ and 4> s ,t- are purely 
imaginary on (a, 6) and satisfy, 

-2<f) sA +(x) = 2<f) sA -(x) = 2m J dv 8i t(u) = g s ,t,+ (x) - g a>t ,-(x), for x G (a, b), (3.12) 

which means that —2<j) S) t and 2<^ s ^ provide analytic extensions of g s t+ ~ 9s,t- into the upper 
half-plane and lower half-plane, respectively. Further, 2g Sy t + 2^> s ,t — V s ,t ~ £s,t is analytic in 
V \ (-00, b] and satisfies by (|3~T2]l and lEOfl) . 

2ff«,t,± + 20 S)t) ± - F S)t - £ S)t = g s ,t,+ + g s ,t- - V S; t ~ t a ,t = 0, on (a, b), 
so that by the identity theorem, 

2g.,t-V ait -£,,t = -2<l>.,t, onV\(-oo,a]. (3.13) 
Using (|3.8|) this yields, 

9s,t,+ + g s ,t- - v s ,t - £ s ,t = 2g s ,t- - V Sjt - £ s ,t + {g s ,t,+ ~ &>,*,-) 

= — 2<f) s ^t- + 2iri, on (— 00, a). (3.14) 

Inserting (|3.12|) . (|3.13|) . and (|3.14f) into (|3.1U() . the jump matrix for T can be written in terms 
of <j) st as 



vt 



v e 2n ^ 
,0 1 



s,t, 



on (a, b), 
on K\ (a, 6). 



(3.15) 



It is straightforward to check, using the fact that <t> a ,t,+ + = on (a, 6), that t;^ has on 
the interval (a, b) the following factorization, 

wr=( e 2ni, i ,_ 1) (_°i 0) ( e 2nl, tl+ J) > on (a, 6), (3.16) 

and the opening of the lens is based on this factorization. Observe that, since Re (f) s ,t,±(x) = and 
Im0 Sitj ±(x) = =F L dv s j(u) for x G (a, 6) (see (|3.12(l ). and since z/ S) t is positive on (a + 5,b — 5) for 
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Figure 1: The lens S 



5 > and s, t sufficiently small (see Remark 12. 3[1 . it follows (as in jSj) from the Cauchy-Riemann 
conditions that 



Re0 Sj t( z ) < 0, for |Im,z| ^ small and a + 5 < Rez < b — 5. 



(3.17) 



We deform the RH problem for T into a RH problem for S by opening a lens as shown in 
Figure so that we obtain a contour S. For now, we choose the lens to be contained in V, but 
we will specify later how we choose the lens exactly. Let 



S(z) 



T(z) 
T(z] 



1 0' 

_ e 2ru/> s ,t(z) i 



for z outside the lens. 

for z in the upper part of the lens, 



(3.18) 



T(z) [„,,.. ) , for z in the lower part of the lens. 
v ' \ e 2n<t> s ,t(z) j / 



Then, using (|3.1(ij) and the RH conditions for T, one can check that S is the unique solution of 
the following RH problem. 

RH problem for S: 

(a) S : C \ S -» C 2x2 is analytic. 

(b) S+(z) = S-(z)vs(z) for z G S, with 



vs = < 



e 2n4> a ,t l J ' 
'-^ e -2n<t> St t,- S 

.0 1 



on (a, 6), 
on S n C±, 
, onR\(a,6). 



(3.19) 



(c) S(z)=I + 0(l/z), 



as z — > oo. 



Remark 3.2 On the lips of the lens (away from a and 6) and on R \ [a — 5, b + 5] , it follows from 
(|3.17|) and (|2.2() that the jump matrix for S converges exponentially fast to the identity matrix 
as n — > oo. This convergence is uniform as long as we stay away from small disks surrounding 
the endpoints a and b. Near these endpoints we have to construct local parametrices. 
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3.4 Parametrix p(°°) for the outside region 

From Remark 13.21 we expect that the leading order asymptotics of Y will be determined by a 
solution p(°°) of the following RH problem. 

RH problem for p(°°): 

(a) P(«0 : C \ [a, 6] -> C 2x2 is analytic. 

(b) Pi 005 (a;) = P^x) f °, J), for x G (a, 6). 



v 1 °y 

(c) P(°°)(ar) = J + C(l/z), as z ^ oo. 
It is well known, see for example (Hl^H; that p(°°) given by 

P^\z) = () l )i( Z -^T lA () Vf\ forzeC\[a,b), (3.20) 



i —i J \z — a t 

is a solution to the above RH problem. Note that p(°°) is independent of the parameters s, t 
and n. 

3.5 Parametrix p(°) near the regular endpoint a 

Here, we do the local analysis near the regular endpoint a. Let = {z G C : \z — a\ < 5} be a 
small disk with center o and radius 5 > sufficiently small such that the disk lies in V. We seek 
a 2 x 2 matrix valued function P^ (depending on the parameters n, s, and t) in the disk Ug, a 
with the same jumps as S and which matches with p(°°) on the boundary a of the disk. We 
thus seek a 2 x 2 matrix valued function that satisfies the following RH problem. 

RH problem for p( a ); 

(a) PW : U s , a \ S -> C 2x2 is analytic. 

(b) P| a) (z) = Pl a) (z)w 5 (z) for z G S n C/ 5 , a , where v 5 is given by lETTO]) . 

(c) p( a ) satisfies the matching condition 

p(«0(js)(p(°°) )(-!)(*) = J + ©(n^ 7 ), (3.21) 
asn-»oo and s, t — > such that ()1.32j) holds, uniformly for z G dUg ia \ £. 

3.5.1 Airy model RH problem 

We will construct P^ by introducing an auxiliary 2x2 matrix valued function _A(£; r) with 
jumps (in the variable () on an oriented contour T = [JjTj, shown in Figure [21 consisting of 
four straight rays 

Ti:argC = 0, r 2 :argC = — , r 3 :argC = vr, T 4 :argC = — — . 

These four rays divide the complex plane into four regions I , II , III , and IV , also shown in 
Figure El Put 

y, = yj (C; r) = Jki {J{C + r)), j = 0, 1, 2, 
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Figure 2: The oriented contour T. The four straight rays Fi, . . . ,T4 divide the complex plane 
into four regions I , II , III and IV . 



with u> = e 3 and with Ai the Airy function, and let, 



A((; r) = V2ne 4 x < 



2/0 -2/2 \ 
y y'o —2/2/ 

-yi -2/2 
-y'\ -v'-ij 

-2/2 yi\ 
-y'2 y'J 

2/0 2/1 \ 

l Wo y'J ' 



for C e I , 
for C G II , 
for C G HI , 
for C € IV . 



(3.22) 



With y'j we mean the derivative of yj with respect to It is well-known, see e.g. |E1^B' that A 
satisfies the following RH problem. 

RH problem for A: 

(a) A is analytic for ( G C\T and for r in C. 

(b) A satisfies the following jump relations on T, 

A + (C) = ^-(C)(_° 1 J), 
A + (C)=A_(C)(J J), 

A+(C)=^-(C)(J J), 

(c) A has the following asymptotic behavior at infinity, 



for c G r 3 , 

for ( G T 1; 

for c e r 2 ur 4 . 



(3.23) 
(3.24) 
(3.25) 



A(C; r) = (C + r)-^iV / + O ( (C + r) 



n-3/2 



e _| (c+r) 3/ 2(J3 



/ - ±r 2 r ^ + ^c 1 / + o(r 6 c~ 3/2 ) + ok- 1 ; 



x e -(§C 3/2 +< 1/2 )- 3 



(3.26) 
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as £ — ► oo, uniformly for r such that 

sgn(lm(C + r)) = sgn(lmC), and \r\ < (3.27) 
In (|3.26|) . N is given by 

N=^=( l , ])e-^ iai . (3.28) 



v^V" 1 1 

3.5.2 Construction of P^ 

We seek P^ in the following form 

pW(z) = E^(z)a 3 A (n 2 ^f a (z);n 2 / 3 r Stt (z^j oge^W" 8 , (3.29) 

where E^ is an invertible 2x2 matrix valued function analytic on U$ a and where f a and r Sj t 
are (scalar) analytic functions on U$ >a which are real on (a — 5, a + 5). In addition we take f a 
to be a conformal map from Us^ a onto a convex neighborhood f a (Us t a) of such that / a (a) = 
and f a (a) < 0. 

If those conditions are all satisfied, and if we open the lens £ (recall that the lens was not 
yet fully specified) such that 

/ a (£n(*7 5 , a nc+)) =r A n f a {Us, a ), and / a (£n (*7 5 , a nc_)) = r 2 n / a (Z7 5 , a ), 

then it is straightforward to verify, using (ETUI]) and lj3~23V(j3~231) . that P<°) defined by (l3~2H|) 
satisfies conditions (a) and (b) of the RH problem for P^ . 



Let 

/«(*) 



2/3 

(a-z) 



as z — » a, (3.30) 



/I , \2/3 

= - U/io(a)V6^J (z-a) + 0((z-a) 2 ), 

where we have used (|2.5|1 . and let 

r s , t (z) = (-mJ\sMt)+tMZ))'%) fa{z)- l/2 . (3.31) 

Then, / a is analytic with f a (a) = and f' a {a) < 0, it is real on (a — 5, a + 5), and it is a 
conformal mapping on U$ a provided 5 > is sufficiently small. Further, it is straightforward to 
check that r s j is analytic on Ug a and real on (a — 5, a + 5), as well. Thus, f a and r S) t satisfy the 
above conditions, so that P^ defined by (|3.29j) . with E^ any invertible analytic matrix valued 
function, satisfies conditions (a) and (b) of the RH problem for P^ . 

Remark 3.3 We can use any functions f a and r Sj j, satisfying the conditions stated under equa- 
tion (|3.29|) . to construct the parametrix P^ a \ However, we have to choose them so as to com- 
pensate for the factor e n<fc > t<T3 in Using (EHTl) . (|3~TB . and the fact that J* ij> 8>t+ (u)du = 1 
we have 

e -n(|/ a ( Z )3/2+r s , t (z)/ a (z)V2) (T3 = (_ 1 )n e -n^ )t ( Z )<73 ) for z £ jj^ \ [ a> a + § ). (3.32) 
From this and 1)3.26 J) it is clear that our choice of f a and r s ,t will do the job. 
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It now remains to determine E^ such that the matching condition (c) holds as well. In 
order to do this we make use of the following result. 

Proposition 3.4 Let n — > oo and s,t — ► such that (jl.32|) holds. Then, 

(n 4 / 7 r st (z)) 2 1/7 (n 4 / 7 r st (z)) 4 2/7 m , „ /7 J , 

•- 4/. ( :;u ^-" /,+ 32/iw j "- 2/7 +o(-- 3/7 )j • <"3) 

Proof. We will use the asymptotics (|3.26j) of A. In order to do this we have to check that 
condition (|3.27f) is satisfied for our choice of £ = n 2 / 3 / a (z) and r = n 2 / 3 r s j(z). 

Obviously r s ^(z) = 0(n~ 4 / 7 ) as n — ► oo and s,t — > such that (|1.32|) holds, uniformly for 
z £ dUs t a- Then, it is straightforward to check that there exists tiq 6 N sufficiently large, and 
K\ , K2 > sufficiently small, such that 

\n 2 '\ S)t {z)\ < \n 2 /*f a {z)\ l '\ (3.34) 

for z 6 dU$ a (for a possible smaller 5), for n > no, and for s and i such that |cin 6 / 7 s — sq| < K i 
and |cin 4 / 7 i — to I < K 2- 

Further, since f a and r s t are analytic near a and real valued on (a — 6, a + 5) one can check 
that 

Im/ a (z) = /^(Rez)Imz + 0((lmz) 2 ), as z — ► a, 



Imr s ^(z) =Tj t (Rez) + C((Imz) ), as z 



a. 



Now, since /J (a) / and r^(Rez) = C(n" 4 / 7 ) uniformly for z e [ft, a one then can find a 
constant C > such that, 

|Imr Sjt (z)| < C|Imz| < |Im/ a (z)|, 

for z G c?f7,5 )a \ (a — 5, a + (5), for n > no, and for s and £ such that |cin 6 / 7 s — so I — K i an d 
|cin 4//7 i — to | < K 2 (for a possible smaller 5, K\ and and for a possible larger no). This yields 

sgn(Im(/ (z) + r Sjt (z))) = sgn(Im/ a (z)). (3.35) 

We now have shown that condition (|3.27|) is satisfied so that we can use the asymptotic 
behavior mj) of A. Using (EOHl) . (l3~2lfl) . (l3~32l . and the fact that r S)t (z) = C(n~ 4 / 7 ) we obtain 
(13331 . □ 

From (|3.33|) and the fact that r S £ = 0(n~ 4 / 7 ) it is clear that (in order that the matching 
condition (c) is satisfied) we have to define E^ by, 

E (a) = (_i)"pM CT3 iV-V3(n 2 / 3 /a) <T3/4 - (3-36) 

Obviously, E^ is well-defined and analytic in Us >a \(a,a + 5). Further, using condition (b) of 
the RH problem for p(°°), equation (|3,28|) . and the fact that = ifa + on ( a > a + <^)> h is 
straightforward to check that -E 1 ^ has no jump on (a, a + 5). We then have that E^ is analytic 
in C/^a except for a possible isolated singularity at a. However, E^ has at most a square root 
singularity at a and hence it has to be a removable singularity. Further, since det P(°°) = 1 and 
det iV = 1 it is clear that det E^ = 1 and thus E^ is invertible. This ends the construction of 
the parametrix near the regular endpoint. 
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3.6 Parametrix near the critical endpoint b 

Here, we do the local analysis near the critical endpoint b. Let Us : b = {z £ C : \z — b\ < 5} be 
a small disk with center b and radius 5 > sufficiently small such U$ t b lies in V and such that 
the disks Ug a and Us b do not intersect. We seek a 2 x 2 matrix valued function P^ (depending 
on n, s and t) in the disk Ug b with the same jumps as S and with matches with on the 

boundary dUgb of the disk. We thus seek a 2 x 2 matrix valued function that satisfies the 
following RH problem. 

RH problem for P^: 

(a) P^ : U s>b \ S -> C 2x2 is analytic. 

(b) P| 6) (z) = P W (2)^(2:) for z G Efyj, n S, where v 5 is given by (EOHl) . 

(c) PW satisfies the matching condition 

p( 6 )( z )(p(°°))- 1 ( z ) = / + 0(n _1/7 ), (3.37) 

as n — > oo and s, t — > such that 1)1.32(1 holds, uniformly for z G dUs )b \ S. 

Due to the singular behavior of the equilibrium measure di>o(x) near 6, see Assumptions II .11 (ii). 
the Airy parametrix does not fit near b. Instead we use a different model RH problem associated 
with the Pf equation (|1.14j) . 

3.6.1 Model RH problem for the Pf equation 

We construct P^ by introducing the following model RH problem for the special solution y of 
the Pf equation (|1.14|) as discussed in Section ll.Ml This RH problem depends on two complex 
parameters s,t and has jumps on the oriented contour T as defined in Section f3. 51 see Figured 
We seek a 2 x 2 matrix valued function = ^(Ci s >t) satisfying the following conditions. 

RH problem for 

(a) * is analytic for £ G C \ T. 

(b) VP satisfies the following jump relations on T, 

* + (C) = MC)(_° 1 J), 
*+(0 = *-(o(J }), 
t + (0 = MC)(I ?)' 

(c) has the following behavior at infinity, 

*(0 = C&N (i - ha 3 C 1/2 + \ C 1 + 0(C 3/2 )) e- e ^°\ (3.41) 

where y = y(s, t) is the special solution of the Pj equation (|1.14l) as discussed in Section 
II ,3| where §j = — y, where N is given by (|3.28j) . and where 9 is given by (|1.21|) . 



for c e r 3 , (3.38) 
for C € Ti, (3.39) 
for C G T 2 U T 4 . (3.40) 
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Remark 3.5 Note that the only difference between the model RH problem for Airy functions 
and the one for P 2 lies in the asymptotic condition (c). In particular, in 6 we have an extra 
factor C 7 / 2 . 

If we fix so,io £ ^ it was proven in [3 Lemma 2.3 and Proposition 2.5] that there exists a 
neighborhood U of sq and a neighborhood W of to such that the RH problem for ^ is (uniquely) 
solvable for all (s, t) G IA x W. Furthermore, for (s, t) £ IA x W, ^ is analytic both in s and t, 
and condition (c) holds uniformly for (s,t) in compact subsets of IA x W. 

In [3 Section 2.3], the authors have shown that the solution ^ of the RH problem for \I/ 
satisfies the Lax pair (fl~ll^ - (fl~Tg|) . From (|L"2"U|) . (|3~1T|) . and (|3~1U)) we then obtain 



$i(C;M) 
*2(C;M) 



*u(C;s,t) 
*2i(C;*,*), 
^n(C;s,t) N 
*ai(C;M), 




for < ArgC < W 7 , 
for 67r/7 < ArgC < 7T - 



(3.42) 



3.6.2 Construction of P( 6 ) 
We seek p( b ) in the following form 

P( b \z) = (n 2/7 / b (z);n 6 / 7 s/ 1 (z),n 4 / 7 t/ 2 (z)) e "*-.*(*)°s, 



(3.43) 



where P^ is an invertible 2x2 matrix valued function analytic on U$b and where /&, /i, and 
/2 are (scalar) analytic functions on Usb which are real on [b — 5,b + 5). We take f\ and f 2 to 
be such that fi(b) = c\ and /2(b) = c 2 (where ci and C2 are given by p.3Ujl ). Then it is clear 
from (|1.32|) that for n sufficiently large and s and t sufficiently small, 



6/7 



sfi(z) G U, and n 4/7 tf 2 (z) G W, for 2: G f7<5 >6 , 



where W and W are the neighborhoods of sq and to where \1/ exists. In addition we take /& to 
be a conformal map from Usb onto a convex neighborhood fb(Usb) of such that /&(&) = and 
/^(fe) > 0. 

If those conditions are all satisfied, and if we open the lens £ (recall that the lens was not 
yet fully specified near b) such that 



/ 6 (En(c/ (5)6 nc + )) = r 2 n/ 6 (c/ 5j6 ), 



and 



/ 6 (sn (E/j 6 nc_)) = r 4 n fb{u s>a ), 



then it is straightforward to verify, using (I3~P31) and (f3~3g|) - (f3~47I| ). that P( 6 ) defined by (^33) 
satisfies conditions (a) and (b) of the RH problem for P^ . 



Let 



AW 



105 



as z — ► 0, where 
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2/7 



,-7/2 



2/7 



6) 



6) + 0(z - bf 



(3.44) 



y/io(6)\/6-a 



To get the expansion of near b we have used (|2.5|) and the facts that /to(&) = ^o(^) = ^ ( see 
(|2.7j) ). Further since /i '(&) > we have that c > 0. So, we have defined an analytic function 



21 



fb with fb(b) = and f' b {b) = c > 0, which is real on (6 — 5, b + S), and which is a conformal 
mapping on Us t b provided 5 > is sufficiently small. 
Next, let f\ and / 2 be defined by 

h{z) = (-m J h Mz)' 1/2 , f 2 (z) = -3 (-mj\(Z)<%) M*)'*' 2 - (3-45) 

Since fb is a conformal mapping in Us,b it is clear from (|2.14|) and Q2.6JI that /i is analytic in 
Us t b- To see that f 2 is analytic in C/^b as well, we also need to use the extra condition h<2.{b) = 
(see (|2.12(l ). Further, f\ and f 2 are real on (b — 5, b + 5) and one can check that, 

Thus, fb, fi, and /2 satisfy the above conditions, so that defined by (|3.43[) . with any 
invertible analytic matrix valued function, satisfies conditions (a) and (b) of the RH problem 
for p0>). 

Remark 3.6 As in Eemark 13.31 we note that we could have also used different functions fb, fi, 
and f 2 . However, we have to choose them so as to compensate for the factor e n<?is ' t<T3 in (|3.43[) . 
Using (fmi . (j33U), $2$ , (E2U), and (|3~TTI) we have 

din 2 ' 7 f b {z)-n & l 7 sh{z),n A / 7 tf 2 {z)) = n^ t (z), for z E \ (6 - 5,6]. (3.47) 

From this and (|3.41|) it is clear that our choice of fb, fi, and fi will do the job. 

It now remains to determine E^ b ' such that the matching condition (c) holds as well. In 
order to do this we make use of the following proposition (the analogon of Proposition 13.4)1 . 

Proposition 3.7 Let n — > oo and s,t — > such that (jl.32|) holds. Then, 

-<73/4 



(3.48) 



where we have used for brevity the notation 

h = h(n 6 / 7 sf 1 (z),n 4 / 7 tf 2 (z)), and y = y(n & ' 7 sh{z),n A l 7 tf 2 (z)). 

Proof. This follows easily from (|3~l3l . (j3~ill . and (|3~47|) □ 

From (J3.48|) it is clear that (in order that the matching condition (c) is satisfied) we have to 
define E^ by, 

E (b) =p (oc) N -lf 2P h Y 3l \ (3.49) 



where ./V is given by Q3.28JI and where p(°°) is the parametrix for the outside region, given by 
(|3.20j) . Similarly as we have proven that £ (a) is an invertible analytic matrix valued function in 
Us, a , we can check that £W is invertible and analytic in U§b- This completes the construction 
of the parametrix near the singular endpoint. 
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Figure 3: The contour Er after the third and final transformation. 



3.7 Final transformation: S i— ► R 

Having the parametrix p(°°) for the outside region and the parametrices P^ and P^ near the 
endpoints a and b, we have all the ingredients to perform the final transformation of the RH 
problem. Define 



R(z) 



' S{z) (PW)" 1 (z), for z G Us, a \ E, 
^(P^)" 1 ^), for z G \ S, 
^S(z) (PM)- 1 (z), for z G C \ (E U U s , a U ^,&). 



(3.50) 



Then, by construction of the parametrices, P has only jumps on the reduced system of contours 
Er shown in Figure El and R satisfies the following RH problem. The circles around a and b 
are oriented clockwise. 



RH problem for R: 

(a) R : C \ Sr -> C 2x2 is analytic. 

(b) P+(z) = R-(z)vr(z) for z G Er, with 



pW(pH)-\ on dU Sia , 
^'(PH)" 1 , on 8Us !b , 
k P (oo) us (PH) _1 , on the rest of Er. 



(3.51) 



(c) P(z) = 7 + 0(l/z), as z -» oo. 

(d) P remains bounded near the intersection points of Sr. 



As n — > oo and s, i — > such that (|1.32j) holds, we have by construction of the parametrices 
that the jump matrix for R is close to the identity matrix, both in I? and L°°- sense on Er, 



^.( z ) 



j + o(„-V7) ) on ac/ 5 , a u ac/ 5i6 , 

I + C(e-T n ), on the rest of Er, 



(3.52) 



with 7 > some fixed constant. Then, arguments as in jlUl lllj guarantee that R itself is close 
to the identity matrix, 



R(z) = I + ©(n" 1 / 7 ), uniformly for z G C \ Er, 



(3.53) 



as n — > oo and s,t — ► such that ()1.32j) holds. This completes the Deift/Zhou steepest descent 
analysis. 
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Remark 3.8 The Deift/Zhou steepest descent method can be generalized to the case where 
the support of v§ consists of more than one interval. However, there are two (technical) differ- 
ences. First, in the multi-interval case, the equilibrium measures v\ and v 2 have densities which 
are more complicated than in the one-interval case, but it remains possible to give explicit for- 
mulae. Consequently, condition (|1.27|) . which expresses the requirement that the density of v 2 
vanishes at the singular endpoint, has to be modified. Further, the construction of the outside 
parametrix p(°°) is more complicated, since it uses 0-functions as in (101 Lemma 4.3]. With 
these modifications the asymptotic analysis can be carried through in the multi-interval case. 

3.8 Asymptotics of R 

For the purpose of proving the universality result for the kernel Kn ( Theorem 1 1.7[) it is enough 
to unfold the series of transformations FwTwSwii and to use (|3.53j) . This will be done in 
the next section. However, in order to determine the asymptotics of the recurrence coefficients 
a*n ,s ^ and 6^ n,s '^ fTheorem ll.ll|) we need to expand the 0(n -1 / 7 ) term in ()3.53j) . 

We show that the jump matrix vr for R has an expansion of the form, 

VR{ Z )=I+^ + ^ + 0(n-^ (3.54) 

as n —* oo and s, t — > such that (|1.32j) holds, uniformly for z G S/j, and we will explicitly 
determine Ai and A2. On S/j \ {dU$ A U 817$^), the jump matrix is the identity matrix plus an 
exponentially small term, so that 

Ai (z) = 0, A 2 (z) = 0, for z G Tir \ (dU s , a U dU s , b ) . (3.55) 

Now, from (I33H) . (l3~3"H|) . and (ETi8l we obtain (1331) with, 

A 1 (z) = -I (n^^^z)) 2 f a {z)- l l 2 P^\z)a^\z)-\ for z G dU s , a , (3.56) 

A 1 (z) = -hf (z)- 1 / 2 P^\z)a 3 P^\zy\ for z G dUs^, (3.57) 



and 



A 2 (z) = i(n 4 \ ( (z)) 4 / s W" 1 /, for z £ dU s , a , (3.58) 

A 2 (z) = lf (z)- 1 ( h2 l y 2 ), for z G dU$ t b, (3.59) 



2 V— iy h 

where we have used for brevity 

h = h{n Q l 7 sh{z)^l 7 th{z)), y = y(n 6 / 7 sf 1 (z),n 4 / 7 tf 2 (z)). 

Observe that Ai and A2 have an extension to an analytic function in a punctured neighborhood 
of a and a punctured neighborhood of b with simple poles at a and b. 

As in 11, Theorem 7.10] we obtain from (|3.54|) that R satisfies, 

m = I+ ^ + ^l + 0(n -^, (3.60) 

as n — > 00 and s, t — > such that ljl.32j) holds, which is valid uniformly for z G C \ (dUs, a U dUs^) ■ 
We have that 

RW and i?( 2 ) are analytic on C \ (dU s>a U dU s>b ), (3.61) 

as z — > 00. (3.62) 
We will now compute the functions R^ and R^ explicitly. 
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Determination of R^ 1 ': 



Expanding the jump relation R + = R-Vr using (|3.54l) and (|3.6U|) . and collecting the terms with 
n -1 / 7 we find 



RX'{z) = R ( _>(z) + Ax(z), for z G dU s , a U dU s , 



b- 



This together with (|3.61j) and 1)3.62(1 gives an additive RH problem for R^\ Recall that Ai is 
analytic in a neighborhood of z = a and z = b except for simple poles at a and b. So, 



Ai(*) 



+ 0(1), as z -> a, Ai(z) 



+ 0(1) 



z — a "2 — 6 

for certain matrices and l^ 1 ). We then see by inspection that 



as z — > 6, 



i?«(z) 



r aC 1 ) 




5« 




+ 




z — a 




z- b 


Ad) 




BW 




+ 




> z — a 




z-b 



ior zeC\(Us, a UUs, b ), 
Ai(z), for z G U s a U U Sb , 



(3.63) 



solves the additive RH problem for Ri 1 ). In now remains to determine A^ and B^K This can 
be done by expanding the formulas (|3.56|) and l|3.57() near z = a and z = b, respectively. We 
then find after a straightforward calculation (using also the fact that /i(6) = c\ and /2(b) = c 2 , 
see 



A« = ^V6^(n 4 / 7 r s ,,(a)) 2 (-/:(a))-^ 



1 j 
i -1 



(3.64) 
(3.65) 



B^ = \hVb^-af b {b)-^^ 1 
where we used h to denote /i(cira 6//7 ,s, c 2 n 4 / 7 i) for brevity. 
Determination of R^: 

Next, expanding the jump relation R + = R-Vr using (|3.54f) and (|3.6U|) . and collecting the terms 
with n -2 / 7 we find 



Rf{z) = R^'(z) + R w (z)A 1 (z) + A 2 (z), for z G dU s , a U dU SJ> . 

This together with (|3.61|) and l|3.62() gives an additive RH problem for R( 2 K Since R W is the 
boundary value of the restriction of Rp-> to the disks Ug >a and Ug^ and since Ai and A2 are 
analytic in a neighborhood of a and b, except for simple poles at a and 6, we have 



>(2), 



i?«(z)A 1 (z) + A 2 (z) 
i?«(z)A 1 (z) + A 2 (z) 



+ 0(1) 



as z — > a, 



z — a 



z-b 



+ C(1), asz 



for certain matrices A^ and B^ 2 \ As in the determination of i?*- 1 ) we then see by inspection 
that 



rM(z) 



( A( 2 ) 


fl(2) 


h 




z — a 


z-6 






h 

v z — a 


z-6 



for z eC\(^U[/ J)t ), 
i?( 1 )(z)A 1 (z) - A 2 (z), for z G Lfr, a U U* b , 



(3.66) 
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solves the additive RH problem for . The determination of and is more complicated 
than the determination of A^ 1 ' and B^ . It involves R^(a) and R^(b) for which we need to 
determine also the constant terms in the expansions of Ai near z = a and z = b. After a 
straightforward (but rather long calculation) we find, 

A m _ i^TsMf ( A {n^r s , t {a)fh fl A 

32(-/i(a)) V"* V 8(-/^(a))V2^( 6) i/2 ^-i lj ' 1 ' 1 

y + h 2 f0 A (n 4 / 7 r s , t (a)) 2 /i /-l 



R ( 2 ) = ^ 1 " / u M V t /t / -1 -t \ 

" 2/^(6) W V Z{-f' a {a)Y/ 2 f' b {b) l/2 V"< -V ' 1 j 

where we used /i and y to denote h(cin 6 ^ 7 s, C2n A / 7 t) and y(cin 6//7 s, C2n 4 / 7 i) for brevity. 



4 Universality of the double scaling limit 



Here, we will prove the univarsality result for the 2-point correlation kernel K~n ' • We do this by 
using the expression (|3.4j) for Kn'^ in terms of Y and by unfolding the series of transformations 

Proof of Theorem II. Ill From equations (|3.4[) . (J3.9|) . and ()3.13j) . the reader can verify that 
the 2-point kernel Kn'^ can be written as, cf. 0E], 

1 



if( s '*)(x,y) = e -n<t>s,t,+ {x) e -n<t> a>t , + 



(y). 



From (j3~l"8|) and the fact that S+ = RP±> on (6 - <5, b + 6) , see (f3~50|) . we have 



27ri(x — y) 
(6) 



ljT+^yjT+Cx) ( Q ) , forx.yeR. 



,(6) 



RPr, on (b,b + 5), 

' 1 0^ 

e 2n,(j> s>t>+ j 



PP 



(6) 



, on (6 — 5, 6). 



Inserting this in the previous equation for iQ ' ; we arrive at, 

= e -^..t.+(*) e -^..*.+ (i/) 1 (0 1) p- 1 (y)R- 1 (y)R(x)P(x) 

lm(x — y) \ u 



for x G (6 — 6, b + 5), where 



r P ( fe ) 



p = < 



,(6) 



1 > 

e 2n(f> S: t,+ y 



on (6, 6 + 5), 
, on (6 — <5, 6). 



Further, we define 



* = < 



1 



on 



on 



(4.1) 



(4.2) 



(4.3) 



where 'F is the solution of the RH problem for ^, see Section 13.61 By (|3.42|) . we have that 
$n = $i and $21 = $2- Using (|3.43j) . (|4.2|) . and (|43|) a straightforward calculation yields, 



P(x) = E^(x)$ (n 2 ' 7 f b {x);n & l 7 sh{x),n A / 7 th{x)) e^'.+ W 



for x G (6 - <5, 6 + 5). 
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Inserting this into (|4.1|) we then obtain, 

K^Hx^) = —± — , (o l) (n 2 /76(y);n 6/ Vi(y),n 4 /V 2 (y)) 

2tti(x — y) \ J 

x (E^)- 1 (y)R- 1 (y)R(x)E^(x)^ (n 2 / 7 ' f b (x);n^ ' S f l {x),n A l 7 tf 2 {y)) Q , (4.4) 

for x £ (b-8,b + S). 

Now, we introduce for the sake of brevity some notation. Let 



u n = b+— m , and Vn = b+^-, with c = f b {b) = (f h^b)y/b=af' ' . (4.5) 

Girt ' cTii ' 

We then have, 

lim n 2 ! 7 fb(u n ) = u, and lim n 2//7 fb{v n ) = v. (4-6) 

n^oo ' n— >oo 

Furthermore, since f\(b) = c\ and /2(b) = c 2 (see 1)3.46(1 ) we have in the limit as n — > 00 and 



s, t -> such that (|L32j) holds, 

limn 6/7 s/i(u n ) = s , limn 6/7 s/i(u n ) = s , (4.7) 

limn 4 / 7 t/ 2 (u n ) = t , hmn 4 / 7 i/ 2 K) = t . (4.8) 

Now, a similar argument as in |24j shows that 

]h^E- l {v n )R{v n )- X R(u n )E b {u n ) = /. (4.9) 

Inserting l(43 |) -((4~U |) into (|4~i)) and using the fact that = <&i and ^21 = ^2 it is then 
straightforward to obtain 

lim —j-K^ (u n ,v n ) 

= ir-7^ — r(o i)tt-Ht>;*o,to)*(u;*o,to)(i 

lni{u — v) \(J 

1 '$i(u;s ,t )^2(t';so,io) - $i(v;s ,to)^2(n;s ,to)) , (4.10) 



2m(u — v) 



where we take the limit n —* 00 and s,i — > such that 1|1.32|) holds. This completes the proof 
of Theorem 11.71 □ 

5 Asymptotics of the recurrence coefficients 

We will now determine the asymptotics of a^' s,t ^ and b^' s '^ as n — > 00 and s,t — > such that 
(|1.32j) holds. In order to do this, we make use of the following result, see e.g. 0^^. Let Y be 
the unique solution of the RH problem for Y . There exist 2x2 constant (independent of z but 
depending on n, s and t) matrices Y\ and Y2 such that 

n*)( Z 7 1) = I+- + ^ + 0(l/z 3 ), asz^oo, (5.1) 



x z n ) z z 2 

and 



it s ' t] = V(*i)i 2 (yi)2i, bt s ' t} = (yi)u + 2^. (5.2) 

{ Y l)l2 
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We need to determine the constant matrices Y\ and Yo.- F° r large \z\ it follows from (|3.9|) . (|3.18|) . 
and (I33U1) . that 



Y{z) = e^^ ta - i R{z)P {oo \z)e n9a ' t{ - z)u3 e-^ nls > 



(5.3) 



So, in order to determine Y\ and Y2 we need the asymptotic behavior of P^°°\z), e n9s ' t ^°" i , and 
R(z) as z — > 00. 

Asymptotic behavior of p(°°\z) as 2 — > 00: 

Expanding the factor ((z - b)/(z - a)) CT3 / 4 in (EHUl) at z = 00 it is clear that, 



p(oo) p(oo) 

pM(z) =/ + ^— + + 0(l/z 3 ), asz^oo, 



with 



(00) 



1 
-1 



R 



'ilr u+ } ( * 1 



(5.4) 



(5.5) 



Asymptotic behavior of e n9a ' t ^ CT3 as z-> 00: 
By (|3.5|) we have 



ng S) t(z)cr a I z 



with 



z 7 



n / udv s ^{u) 



G G 

I + — + ^- + 0(l/z 3 ), asz^oo, 



z z £ 



(5.6) 



1 

-1/ ' 



G 2 



* 
* 



(5.7) 



Asymptotic behavior of R(z) as z — > 00: 

As in the matrix valued function R has the following asymptotic behavior at infinity, 

R(z) = I + ^l + ^l + 0(l/z 3 ), asz^oo. 
z z l 

The compatibility with ()3.6(Jj) . ()3.63|) . and (|3.66f) yields that 

Ry = (AW + 5(D) „-V7 + ^(2) + B (2)j n -2/7 + ( n -3/7 }j 
ij 2 = ( a ^(D + bB (D) n -V7 + ( aA (2) + bB {2)\ n -2/7 + (n~ 3 / 7 ), 



(5.8) 
(5.9) 

as n — » 00 and s,i — > such that (|1.32|) holds. Here, AW, Z^ 1 ), A^ 2 ), and -B^ 2 ) are given by 
(ESU), (l3~K5l . (I3~F71) . and (l3~HHl) . respectively. 

Now, we are ready to determine the asymptotics of the recurrence coefficients. 
Proof of Theorem ITTTT1 Note that by (l5~H . (ETfijl and (Oil . 



K = e2 



P^ + Gx + Rt 



e 2 



and 



r 2 = e2 



P 2 (oo) + G 2 + i? 2 + ^ 1J P 1 (oo) + (p^ + i?! ) Gi 



e 2 



kni s ,tcr3 



(5.11) 



(5.12) 
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We start with the recurrence coefficient a£ l ' s,t \ Inserting Q5.11[) into (|5.2|) . and using the 

facts that {P^^n = -{P^hi = i(b-a)/A (by (JS3J) and (G x )i 2 = (Gi) 2 i = (by (52)) , we 
obtain 



7 (n,s,t) 



' j " ^ + ((^1)21 " (#1)12) + (#1)12(^1)21 



1/2 



(5.13) 



Now, from the formula ()5.9j) for Ri and the formulas (|3.64j) . (|3.65() . Q3.67JI . and (|3.68j) for 

AW,bW,A<V, and 

, we have 



(#1)21 — (#1)12 



y / ( n 4 /V st (a)) 2 ^ 



+ 



+ 



#(&) U(-/»)V2 



-2/7 + 0(n ~3/ 7)) 



and 



(#1)12(^1)21 



b - a ( (n 4 / 7 r Sjt (a)) 2 h 



+ 



^-faia)) 1 ' 2 f' b (b) 



n 



-2/7 



+ 0(n- 3 / 7 ). 



Note that we have used y to denote y(cin 6 / 7 s, c 2 n 4 / 7 t) for brevity. Inserting the latter two 
equations into (|5.13|) and using the fact that f^{b) = c (by (|3.44|) 1 we then obtain (|1.39|) . 

We will now consider the recurrence coefficient &i n,s '*' ) . Inserting (|5. 1 1|) and ()5.12j) into (|5.2|) . 

and using the facts that (P-f^n = (P 1 { °° ) )22 = 0, (Gi) w = (G 2 )i2 = 0, (Gi)n + (Gi) 22 = 0, 
and Ri = 0{n~ 1 / 7 ), we obtain 



b (n,s,t) = (R) (# 2 (0O) )l2 + (#l)ll(Pr J )l2 + (#2)l2 

(P^)l2 + (#l)l2 



,(00)-, 



(#l)ll + 



J 12 

l-^l J 12 



+ (#l) 



11 



(A 



2J12 



{Pi J 12 



(#1) 



12 



(H 



1)12 



[Pi J 12 



Am J12, 



+ 0(n- 3 / 7 ) 



Since (P^ia = *(& - o)/4, (Pj°° j )i 2 = i(6 2 - a 2 )/8, i?i = ©(n- 1 /?), and r 2 = 0{ n - l l r ) we 
then obtain after a straightforward calculation and combining terms, 



„(„,,,) = ^ + f 2{fll)ll+2i ^ (iJl)l2 « 

2 V b — a b — a 



(#2)12 
U 



1 + 



4i 



(R 



1 12 



b — a 



b — a 

(Ri)n(Ri)i2 + 0(n-V 7 ). (5.i4) 



Now, from (|5~3J) . (tQ71|) . (jSH, dSH, (3221), and (l3~?^|) we have 



2(^0!! + 2i h -±^{R l ) 12 - -r^(R 2 )i2 = 2i 
b — a b — a 



4(2) _ R (2) 

^12 -°12 



n 



-2/7 



+ 0(n- 3 / 7 ) 



+ 



h 2 (n*/ 7 r s , t (a)) 



f b (b) K(b) 16(-/») 



-2/7 



+ £}(rT 3 / 7 ), 



29 



and 



(i?l)ll(i?l)l2 



A 



B 



n 



-2/7 + 0(n -3/ 7) 



b-a h 2 



A/7 



h? _ (n^ 'r s ,t(a)Y 
J' b {b) 16(-/£(a)) 



n 



-2/7 



+ 0(n- 3 / 7 } 



Inserting the latter two equations into (|5.14j) and using the facts that (i?i)i2 = 0(n and 
/&(&) = c we obtain 1)1.40(1 . So, the theorem is proven. □ 
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